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Brownian vortexes are stochastic machines that use static non-conservative force fields to bias 
random thermal fluctuations into steadily circulatingcurrents 0.i- The archetype for this class 
of systems is a colloidal sphere in an optical tweezer p, ii- Trapped near the focus of a strongly 
converging beam of light, the particle is displaced by random thermal kicks into the nonconservative 
part of the optical force field arising from radiation pressure [^, which then biases its diffusion 
[I|, [^. Assuming the particle remains localized within the trap, its time-averaged trajectory traces 
out a toroidal vortex. Unlike trivial Brownian vortexes, such as the biased Brownian pendulum, 
which circulate preferentially in the direction of the bias, the general Brownian vortex can change 
direction and even topology in response to temperature changes. Here we introduce a theory based 
on a perturbative expansion of the Fokker-Planck equation for weak non-conservative driving. The 
first-order solution takes the form of a modified Boltzmann relation and accounts for the rich phe¬ 
nomenology observed in experiments on micrometer-scale colloidal spheres in optical tweezers. 


I. INTRODUCTION 

Stochastic systems can be driven out of equilibrium in 
various ways. Thermal ratchets and Brownian motors 
use time-dependent forces to rectify Brownian motion 
Other systems evolve in response to spa.tial and 
temporal variations in the temperature [lOl lll| . Even 
if none of these systems can reach thermodynamic equi¬ 
librium, some can achieve non-equilibrium steady states. 
Recently, a distinct class of stochastic machines was dis¬ 
covered that use noise to extract work from a static non¬ 
conservative force field These ^sterns, which 

have been dubbed Brownian vortexes P, 0 , perform no 
work at all in the absence of stochastic forces. When 
activated by noise, they enter into steady-state motion 
characterized by toroidal vortexes in the time-averaged 
or ensemble-averaged probability current. These cyclic 
processes in principle can be coupled to external systems 
to extract useful work. Since their initial observation 
in the fluctuations of optically trapped colloidal beads 
P, [l^ , Brownian vortexes have been reported in trapped 
colloidal rods [l^, in trapped spheres subjected to shear 
flows [a la. i n the response of bacterial populations to 
antibiotics [l6[, and in models for population dynamics 

pa. 

Brownian vortexes arise in time-independent force 
fields that include at least one point of stable equilib¬ 
rium. In the absence of thermal forces, such systems re¬ 
main stationary at their fixed points and so perform no 
work. They differ in this respect from conventional ma¬ 
chines that move deterministically under the influence of 
non-conservative forces. Random thermal forces allow a 
Brownian vortex to explore its force landscape. Were the 
force field purely conservative, the system would reach 
thermal equilibrium in the Boltzmann distribution, and 
so would have no means to perform work. In force fields 
with solenoidal components, however, the probability dis¬ 
tribution can be advected by the non-conservative force. 
Probability currents then flow through the system un¬ 
der the competing influences of advection and diffusion. 


These currents must form closed cycles if the system’s 
overall probability is conserved, raising the possibility 
that the system can reach steady state. 

Initial studies have identified two broad categories of 
Brownian vortex behavior: trivial Brownian vortexes 
that circulate in the direction dictated by the non¬ 
conservative component of the force landscape, and gen¬ 
eral Brownian vortexes that select their own topology 
and circulation. Both types of behavior have been ex¬ 
plained with master equations on discrete networks 0. 
Solutions for continuous systems have been obtained for 
special cases [1, [Tdl - fTsIl , with most examples representing 
trivial Brownian vortexes. Here, we introduce a pertur¬ 
bative theory that accounts for both trivial and general 
cases in continuous systems. 


II. THEORY 

Our approach is based on the observation that systems 
governed by linear forces can be mapped to a general¬ 
ization of the Ornstein-Uhlenbeck process, for which a 
general solution is known [T^ . For nonlinear forces, a 
common approach is to solve the system using perturba¬ 
tion theory [^ . Here we develop a perturbation theory 
for the specific case where the non-conservative force is 
weak compared to the conservative one. We show that 
this approach captures the important features of general 
Brownian vortex circulation, including topological tran¬ 
sitions and flux reversal. 

Following earlier studies P, 0, [H, [13 , we describe a 
system capable of undergoing steady-state stochastic cir¬ 
culation as a Brownian particle of mobility /x moving 
through a static force landscape F(r) at absolute tem¬ 
perature T. Such a description applies naturally to the 
motion of a colloidal particle in an optical force field, 
and also may be applied to more general systems such 
as an ensemble of identical particles interacting through 
conservative forces and confined by a force landscape. 
Assuming that the system reaches steady state, we seek 
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the steady-state probability distribution p(r) describing 
the likelihood of finding the particle near position r, and 
the associated steady-state flux of probability j(r) flowing 
through that point. The flux is generated both by advec- 
tion of the probability distribution and by diffusion: 

j(r) = MpWFW - T>Vp(r). (1) 

Because the probability density is non-negative, we may 
express it as the exponential of an effective potential ^(r): 

p(r)=e-^'^W, (2) 

so that 

j(r) =/ip(r)[F(r)-HV(/>(r)]. (3) 

The Helmholtz decomposition theorem guarantees 

that the force field can be separated into a conservative 
gradient term and a non-conservative solenoidal compo¬ 
nent, 

F(r) =-V17(r)-h V X A(r), (4) 

where U (r) is the potential energy and A(r) is the vector 
potential. The probability current may be written in 
terms of these potentials as 

j(r) = /rp(r) {V x A(r) -t- V [<^(r) - U (r)]} . (5) 

Equation m makes clear that the system can only 
reach thermodynamic equilibrium with jo(r) = 0 if the 
non-conservative force vanishes, V x A(r) = 0. The re¬ 
maining term in Eq. ([5|) then yields Boltzmann’s distri¬ 
bution for the probability density, 

Po{y) = (6) 

with 0o(r) = pU{r), where = ksT is the thermal 
energy scale. 

Systems subject to nonconservative forces may not 
reach thermodynamic equilibrium, but still must satisfy 
the continuity equation. 


whether or not it is conservative. Eor a potential force, 
the Einstein relation D = ksTp, between the diffusion 
coefRcient D and mobility p follows from the condition 
that there exists an equilibrium state which obeys the 
detailed balance condition j(r) = 0. In the presence of a 
non-conservative force, the system does not come to equi¬ 
librium, it does not satisfy detailed balance. Neverthe¬ 
less, its probability distribution remains independent of 
the particle’s transport properties, and depends only on 
the form of the force field. In this sense, the steady-state 
probability distribution resembles Boltzmann’s distribu¬ 
tion despite the system’s departure from equilibrium. 

To find the specific steady-state probability distribu¬ 
tion, we introduce the projection 

p(r) = VC/(r) • V X A(r) (9) 

of the non-conservative force onto the direction of the 
conservative force. An exact solution to Eqs. m and 
(|S]) is known Q only for the special case p(r) = 0. In 
this case, the steady-state probability distribution retains 
the form of the Boltzmann distribution, Eq. (0), and is 
simply advected by the nonconservative force: 

j(r) = X A(r). (10) 

The direction of j(r) in this case is fixed by V x A(r) 
regardless of the temperature. Equation (nni) therefore 
describes a trivial Brownian vortex. Examples of trivial 
Brownian vortexes include the biased Brownian pendu¬ 
lum and colloidal spheres circulating in circularly polar¬ 
ized optical tweezers [2l|. 

To move beyond this special case, we consider systems 
in which V x A(r) is not necessarily aligned with VU (r), 
but may be treated as a perturbation whose scale is char¬ 
acterized by a small parameter e. We therefore replace 
V X A(r) in Eq. (|3]) with eV x A(r). Assuming the per¬ 
turbation not to be singular, the effective potential may 
be expanded in orders of e as 

(j){r) = <j)o{r) + +0{e^} , (II) 


V-j(r) 


dpir,t) 

dt 


with 00 (r) = U{r). The associated expansion of the 
(7) probability current, 


where p{r,t) is the time-dependent probability density. 
Equation © is the Fokker-Planck equation for Brownian 
vortexes. Any steady-state solution of Eq. ([7]) satisfies 

V-j(r)=0, (8) 

which also ensures conservation of probability. 

Quite remarkably, Eq. (|S]) implies that the steady state 
distribution p(r) does not depend on the particle’s mo¬ 
bility, p, or on any other transport property. Indeed, 
substituting Eq. (O for the current density into Eq. (|8]) 
yields an equation that is independent of p. This hap¬ 
pens because of the assumption, implicit in Eq. o, that 
the particle responds in the same way to an applied force 


j(r) = eji(r)+0{e2}, (12) 

has no contribution at zero-th order in e. The first-order 
correction, 

ji(r) = [V X A(r) -h V0i(r)] (I3a) 

retains the advective term from Eq. m, although this 
now distorts the probability distribution as well as trans¬ 
porting it. The second term in Eq. (I13al) describes diffu¬ 
sive relaxation of that distortion. 

Because these two terms depend on temperature in 
different ways the first-order expansion admits the possi¬ 
bility of temperature-dependent transitions such as those 
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(a) (b) (c) 



FIG. 1. (color online) (a) Force field F(r) experienced by a 1.0 jim diameter colloidal silica sphere of refractive index 1.45 in 
a 1 mW Gaussian optical tweezer propagating upward along the z axis with a 60° convergence angle. Golors indicate relative 
intensity of the beam. Arrows denote direction of net force in the (r, z) plane. The superimposed semicircle indicates the size 
of the colloidal particle, (b) Axial component of the restoring force, Fz{z), along r = 0. (c) Radial component of the restoring 
force, Fr{r), in the plane z = 0. Shaded rectangles in (b) and (c) indicate the particle’s range of motion within the force field 
as a function of fluctuation energy in steps of 2 ksT up to 10 fcsT. (d) Non-conservative component of the force, Fz{r) in the 
plane z = 0. (e) Detail of (d) showing fit to a quartic polynomial. 


characterizing general Brownian vortexes. These transi¬ 
tions, moreover, may be understood to arise from com¬ 
petition between advection and diffusion. The balance 
depends on the distortion of the probability distribution 
away from /3o(r). 

As explained in Appendix substituting Eq. (I13al) 
into Eq. m yields a differential equation for (pi (r), 

V^(^i(r) — P'VU{r) ■ V(()i(r) = f3p{v). (13b) 

Solutions to Eq. (I13bl) are difficult to find for arbitrary 
force fields. Provided that C/(r) has at least one min¬ 
imum, however, the associated operator Hjj = — 

/3Vt/(r) • V has eigenfunctions with eigenvalues 

Xn that satisfy 

Huipnir) = Xn'ipn{r). (13c) 

Symmetrized functions of the form consti¬ 

tute a complete orthogonal basis with the orthonormal¬ 
ization condition 

J dr = Snm- (13d) 

We therefore can expand (/>i(r) in this basis, 

(pi{r) = ^ Cn'ipnir) (13e) 

n 


with expansion coefficients 

Cn = ^ /e“^^('’Vn(r)p(r)dr. (13f) 

J 

Equations (TT^ are the principal result of this study. 
They describe the lowest-order extension of a trivial 
Brownian vortex into a general Brownian vortex under 
the influence of a weakly non-conservative force. The na¬ 
ture of the resulting topological transformations and flux 
reversals depends on details of F(r), and particularly on 
the projection of its solenoidal component onto its con¬ 
servative part. 

To illustrate practical applications of Eq. m, we next 
apply this formalism to the particular case of colloidal 
spheres moving in optical tweezers, the context in which 
Brownian vortexes were first observed P, Q . This analy¬ 
sis casts new light on the nature of this system’s behavior. 


III. COLLOIDS IN AN OPTICAL TWEEZER 

To facilitate comparisons with experimental studies of 
colloidal spheres circulating in optical tweezers, we nu¬ 
merically compute the forces acting on a trapped sphere 
with the Lorenz-Mie theory of light scattering using 
methods described in Appendix |B] Typical results are 
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presented in Fig. [TJ Colors in Fig. [ija) represent the rel¬ 
ative intensity of the light in an optical tweezer, as viewed 
in the (r, z) plane in cylindrical coordinates. The light 
propagates in the +z direction (upward) and comes to a 
focus along the axis defined by r = 0. Arrows indicate 
the direction and strength of the resulting force F(r) ex¬ 
perienced by a colloidal sphere of radius Up = 0.5 pm at 
position r within that light field. We have shifted the 
origin of the coordinate system to coincide with the posi¬ 
tion of the trap so that F(0) = 0. Away from the origin, 
F(r) directs the particle back to the stable fixed point. 
All forces reported in Fig. [T] are computed for a power 
of 1 mW, which is a reasonable scale for typical optical 
trapping experiments. 

Both the axial component of the force, plotted in 
Fig. HKb), and the radial component plotted in Fig. [TKc) 
resemble a linear restoring force over a reasonably wide 
range of axial displacements. The axial component of the 
optical force shows a non-trivial dependence on radial po¬ 
sition, as shown in Figs.[T](d) and (e), that is reasonably 
modeled as a quartic polynomial. This solenoidal depen¬ 
dence is responsible for Brownian vortex circulation in 
optically trapped colloidal spheres. 

The three-dimensional force field is very nearly sym¬ 
metric with respect to rotations about the z axis. Small 
distortions along the axis of the light’s polarization have 
little influence on the trapped particle’s motions, and will 
not be considered here. 


A. Simulated circulation 

The data in Fig. Ha) through He) show results of 
Brownian dynamics simulations [22| of a colloidal sil¬ 
ica sphere diffusing through water at T = 300 K in the 
computed force field from Fig. [ija)- The sphere’s tra¬ 
jectory rp(t) is computed at 100 ps intervals for each 
of three values of the light’s intensity: 1.38 mW in 
Fig.Ha), 0.35mW in Fig.^b) and 0.17mW in Fig.^c). 
These correspond to values of the radial stiffness, k, 
of lpNpm“^, 0.25pNpm“^ and 0.12pNpm“^, respec¬ 
tively. These values are chosen for consistency with pre¬ 
vious experimental studies UM- Each trajectory then is 
compiled into estimates for the probability density p{r) 
and the probability current density i(r) using an adap¬ 
tive kernel density estimator [I|, y, ll^. Streamlines of 
j(r) are denoted by barbs whose sharp ends point in the 
direction of motion and whose length is proportional to 
the local current density. Each barb is colored by the 
estimate for p{r). 

At the highest intensity. Fig. Ha); the trap is stiffest 
and the particle is most strongly confined. A single 
toroidal roll is evident in the probability current, and 
circulates downstream along the optical axis. Reducing 
the intensity weakens the trap and frees the particle to 
explore more of the force landscape. Under these con¬ 
ditions, shown in Fig. Hb), streamlines of j(r) reveal a 
second counter-rotating roll. Reference [I| identifies the 
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FIG. 2. (color online) Streamlines of the current density for a 
1 pm-diameter silica sphere in an optical tweezer as a function 
of trap intensity. Sharp ends of the symbols point in the direc¬ 
tion of j(r) with colors determined by the relative probability 
p(r), as indicated by the inset color bar. Light propagates in 
the z direction, (a), (b) and (c) show results from Brownian 
dynamics simulations in force fields computed for a Gaussian 
optical tweezer at a vacuum wavelength of 532 nm and a nu¬ 
merical aperture of 1.4. (d), (e) and (f) show corresponding 
results of analytical expressions in Eq. (EZI for the same pa¬ 
rameters. (a), (d) ^/3kap = 30.5: a single toroidal roll is 
visible near the optical axis, (b), (e) = 7.6: concen¬ 

tric counter-rotating toroidal vortexes, (c), (f) ^pka^ = 3.6: 
weak confinement showing a single flux-reversed roll. All of 
the qualitative features of the Brownian vortex circulation 
observed in the simulation results are obtained also in the 
analytical theory. 


appearance of this second roll as a topological transition. 
Reducing the intensity still further allows the second, 
outer roll to dominate the system’s dynamics. It sub¬ 
sumes the inner roll so that only a single toroidal vortex 
remains. The single remaining roll circulates with the 
flux directed upstream along the optical axis, which sig¬ 
nals a flux reversal relative to Fig. H^-) addition to 
the topological transition relative to Fig. Hb). Precisely 
this behavior was reported in [ij, which lends credence 
both to the experimental results and also to the present 
simulations. 

The particle’s residence time in the trap is effectively 
indefinite at the highest laser intensity. This corresponds 
to the force held in Fig.[T]and the streamlines in Fig.H^)- 
Reducing the laser intensity reduces the residence time, 
which falls to 10s under the conditions in Fig. He)- This 
is consistent with results of experimental studies. Sim¬ 
ulations under these conditions are restarted every time 
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the particle escapes, and the results averaged to until the 
computed current density converges. The streamlines in 
Fig. [SJc) therefore should be viewed as an ensemble av¬ 
erage. 


B. Analytical model 


The force field F(r) resembles a cylindrically symmet¬ 
ric harmonic well for small excursions away from the equi¬ 
librium point. This can be seen in Figs.[ljb) and (c). We 
therefore model the trapping potential as 

U{v) = ^k{r^+r]z^), (14) 

where rj characterizes the trap’s anisotropy. In this case, 
Eq. (I13cl) has a complete set of eigenfunctions, given in 
cylindrical coordinates by 

ipnir) = , (15a) 


where Tn(0 is the Laguerre polynomial of index n, Hn^ (•) 
is the Hermite polynomial of index and where n = 
{n, Tiz} is a set of whole-number indexes that are related 
to the associated eigenvalues by 

Xn =—Pk{2n + rjnz). (15b) 


The harmonic well’s basis functions are normalized by 


Nn = 



1 


(15c) 


We have chosen basis functions that are independent of 
9 to reflect the azimuthal symmetry of F(r). 


1. Quadratic perturbation in an isotropic trap 

To illustrate applications of Eq. (US, we first consider 
the symmetric case, r] = 1, subjected to a quadratic per¬ 
turbation, 

Fz{r) =ekap^l-^^ , (16) 

with Fx{r) = Fy{r) = 0, where e is the small param¬ 
eter characterizing the strength of the non-conservative 
force and Op is the radius of the sphere. This perturba¬ 
tion is divergence-free and thus has the solenoidal form 
assumed in Eq. Choosing Fz{r) to be proportional 
to k ensures that it scales with the light’s intensity in 
the same way as the linear restoring force. Scaling dis¬ 
tances by the particle’s radius is reasonable because typ¬ 
ical realizations of Brownian vortex circulation involve 
motions substantially smaller than Op Hi- In general, 
both k and e depend on the particle’s radius relative to 


the wavelength of light. Parameterized in this way, the 
perturbation may be considered weak if 


e < \I^Pkaj. (17) 

Choosing a quadratic form for Fz{r) suggests that the 
non-conservative force will dominate the linear restoring 
force at large distances, and that the particle eventually 
will escape from the trap. The rate at which probability 
leaks from the system is limited by the smallness of e. 
This is the case, for example, in the system described by 
Fig. [H Under these conditions, the particle is likely to 
remain trapped over the period of observation, and prob¬ 
ability may be considered to be conserved. This also is 
the case for the exTOrimental studies of optically trapped 
colloidal particles |l|,[l,0 that Eqs. m and m are in¬ 
tended to model. 

The projection factor for this system, defined in 
Eq. ([2]), has the form 


p{r) = -k'^GpZ ^1 - . (18) 

It does not explicitly include e because this parameter is 
introduced in Eq. (ITTI) for the effective potential and in 
Eq. (1121) for the current density. Substituting this form 
for p{r) into Eq. (1131) along with the basis functions from 
Eq. (fist yields the first-order correction to the effective 
potential, 


/3<(’i(r) 


3 


4 — Pkttp 



(19) 


This correction’s linear dependence on z displaces the 
probability distribution down the optical axis, as would 
be expected of radiation pressure. The distribution’s 
width 


^p{r)dr 


rp{r)dr 


= \hz ( 20 ) 


2 = 0 


is the same as in the unperturbed case. 
The associated current density, 


j(r) = -eDe ^jdkrzf + (2 — f3kr^)z] , (21) 

o 


has the form of a toroidal roll, centered on the axis r = 0 
and circulating around a circular core in the plane z = 0 
at radius 



which coincides with the distribution’s width, tq = a. 
The circulation rate is controlled by the strength e of the 
driving force and the rate at which the particle diffuses 
in the harmonic potential energy well given its diffusion 
coefficient D = pksT. 

This Brownian vortex undergoes no topological tran¬ 
sitions or flux reversals. Applying a quadratic perturba¬ 
tion to a harmonic well therefore constitutes a model for 
a trivial Brownian vortex. 
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2. Quartic perturbation in an anisotropic trap 

More interesting behavior arises in a more highly struc¬ 
tured force field of the type actually observed in optical 
trapping experiments [if, H, 0- Figure [IJe) shows that 
the computed axial force is well approximated near the 
plane z = 0 by a quartic polynomial, 

F^r) ={l - p^r'^) , (23) 

O 


by introducing a super-Gaussian tail at large r. This 
broadening may be of interest for precision measurements 
of optical forces because it can influence 0 calibration 
protocols based on analysis of thermal fluctuations Hi- 
l26[. The availability of an analytical form for ()>i(r) will 
help to assess when nonequilibrium redistribution of the 
probability density may be ignored (2^ . 

The first-order correction to the current density, 

jiW = jr(r)f-h j^(r)£, (27a) 


with Fx{r) = Fy(r) = 0, where p 4 characterizes the par¬ 
ticle’s interaction with the focused beam of light. Like 
e and fc, this additional parameter also depends on the 
sphere’s radius, ap. As for the quadratic case, Eq. (l23l) 
describes a purely solenoidal perturbation. We will show 
that Pi governs topological transitions and flux reversals 
in the resulting Brownian vortex circulation. Incorporat¬ 
ing a constant offset into Fz(r) would move the equilib¬ 
rium position along the optical axis as in the quadratic 
case, but does not otherwise influence the system’s be¬ 
havior. We have omitted such an offset from Eq. (1^51) for 
clarity. The weak-perturbation condition for this model 
is 


Accounting for the trap’s anisotropy p allows for com¬ 
parison with optical trapping experiments. The force 
field presented in Eig.[Tl for example, is intended to model 
the influence of Gaussian beam brought to a diffraction- 
limited focus and is characterized hy p = 0.47. Results 
for isotropic harmonic wells can be retrieved by setting 
p = 1. No qualitative features of Brownian vortex circu¬ 
lation depend on the value of p. 

The projection factor, defined in Eq. m, associated 
with this force field 

pW = —(^'^ -P4r^), (25) 

O Ojp 

gives rise to a first-order correction to the effective po¬ 
tential that also has the form of a quartic polynomial: 


0 dp 

(1 2pi){2 + r‘^) piT^ 

4:7} 

(26a) 

For conciseness we have introduced the dimensionless pa¬ 
rameters 


= —pj3kr^ and 

(26b) 


{4 + p)pk^^’ 

(26c) 


As in the quadratic case, the quartic perturbation dis¬ 
places probability along z. The quartic term also tends 
to broaden the probability distribution relative to po{v) 


has the radial component 
>(r) 


Jo(i-) 


= pfjkzr 




(27b) 


that is conveniently measured in units of an overall scale 


Jo(r) = 




3 Oh 


(27c) 


that again reflects diffusion in the potential energy well, 
as in Eq. (EU. The radial component of the current 
density changes sign as it passes through the equatorial 
plane, z = Q. The axial component. 


= 2p4^-^f'‘-4 -^^ 4 ) 7^-k4(l-2p4), (27d) 

Jo(r) ry VP / 

varies non-monotonically with distance r from the axis. 
The net current, j(r) = eji(r), is proportional to the 
strength of the non-conservative driving force. 

Streamlines of j(r) plotted in Eig. [JKd) through [JKf) 
reveal circulatory flows that agree well with the simu¬ 
lated results for the same system under the same condi¬ 
tions. Specifically, the first-order result for the strongly 
confined system in Fig. [2Kd) shows a single roll consis¬ 
tent in position, extent and speed with the simulation 
result in Fig. [2Ka). The double-roll structure in Fig. [5Ke) 
similarly is consistent with that in Fig. [2Kb), although 
quantitative details of the circulation differ, particularly 
near the optical axis. The single-roll structure in Fig.j^Kf) 
again qualitatively resembles that in Fig. He), although 
the centers of circulation appear at different radial posi¬ 
tions. The first-order perturbation theory developed in 
Eq. (fOl) thus captures the essential features of general 
Brownian vortex circulation in this model system, albeit 
with quantitative discrepancies. Trends in the analytical 
results therefore offer useful insights into the nature of 
the phenomenon. 


C. Topological transition and flux reversal 

The probability current, j(r), vanishes at the cores of 
toroidal vortexes. Solutions of ji(r) = 0 take the form of 
circles in the plane z = 0 centered on the axis at r = 0. 
Two vortex cores exist if pi is sufficiently small, at radii 
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r_|_ and r_ that satisfy 


r\. = Tg (1 ± A^), where (28a) 

f3k2 + ri \pi ) 

^2 ^ V (8 + 4?7 + 7?2)p| - 4 p 4 +1 
1 + ?7P4 

According to Eq. (l28ll . the vortex at r+ is present for 
all temperatures greater than zero and moves outward as 
the temperature increases. The other vortex at r_ moves 
inward, and ceases to exist when it reaches r_ = 0, which 
occurs when p 4 = i. This condition corresponds to a 
threshold temperature 


ksTc 


4 + r; k 

16 P4 


(29) 


below which the probability current consists of two 
counter-rotating toroidal vortexes and above which only 
a single roll remains. The threshold temperature is pro¬ 
portional to the light’s intensity through the trap stiff¬ 
ness, k, and depends on details of the particle-light in¬ 
teraction through p 4 . Interestingly, it does not depend 
on the strength, e, of the non-conservative force. In most 
optical trapping experiments, the temperature remains 
constant while k is varied. Equation (0^1) for the thresh¬ 
old temperature then may be recast into an equivalent 
condition for the trap stiffness. 

The system’s behavior for T < corresponds to 
strong confinement by the harmonic well. The outer 
toroidal roll therefore may not be perceptible in an ex¬ 
periment of finite duration because the particle spends 
comparatively little time exploring the outermost reaches 
of the force landscape. The apparent transition from one 
vortex in Eig.[2]^a) to two in Eig. [21(b) in fact results from 
increasing the occupation of the outer roll and does not 
constitute a topological transition. The appearance of 
a second concentric vortex in previous optical trapping 
studies Hi almost certainly corresponds to the same 
statistical sampling considerations. 

Increasing the temperature beyond Tc, oi equivalently 
reducing the stiffness of the trapping potential, causes 
an actual topological transition by eliminating the inner 
roll at r_. In this case, the single remaining roll cir¬ 
culates upstream along the optical axis. The behavior 
in Fig. |2](f) therefore reflects a flux reversal relative to 
Fig. |2Id) as well as a topological transition. This more 
dramatic change also has been observed in experimental 
studies [l|. 


D. Size-dependent crossovers 

The nature of the Brownian vortex behavior in an op¬ 
tically trapped colloidal sphere depends qualitatively on 
the form of the nonconservative force, Fz{r). Neither the 
concentric roll structure nor the high-temperature topo¬ 
logical transition appear if this force is quadratic rather 



FIG. 3. (color online) (a) Dependence on the particle radius 
of the quadratic and quartic coefficients e and P 4 character¬ 
izing the non-conservative optical force F^{r). These coeffi¬ 
cients are obtained by polynomial fits to the force field com¬ 
puted with generalized Lorenz-Mie theory for a silica sphere 
in a Gaussian optical tweezer. Results are shown for a silica 
sphere in water and an optical tweezer created with a vacuum 
wavelength of 532 nm in a beam with a convergence angle of 
60°. (b) Fz{r) for a sphere with ap — 0.13pm, together with 
a quartic fit. (c) Cp = 0.45 pm. (d) ap — 0.5 pm. 


than quartic. What dynamic patterns emerge therefore 
depends sensitively on the size and composition of the 
sphere and the structure of the light. 

The data in Fig.|3|show how the controlling parameters 
e and p 4 vary with sphere radius ap for a colloidal silica 
sphere trapped in water by an optical tweezer of vacuum 
wavelength 532 nm brought to a focus by a lens of numer¬ 
ical aperture 1.4. Figures |3Kb), (c) and (d) show plots of 
Fz{r) for the representative radii indicated by vertical 
dashed lines in Fig.|3Ka). As in Fig. [U these results were 
obtained numerically using Eq. (IB6I) . 

In the Rayleigh limit, for ap -C A, the non-conservative 
force is peaked at the optical axis and falls off very nearly 
quadratically in r. For such particles, e > 0. Figure |3Kb) 
is representative of this range of conditions. The quar¬ 
tic contribution being weak, Rayleigh particles enter into 
single-roll Brownian vortexes circulating down the optical 
axis in the direction of the light’s propagation. This is the 
mode of operation that first was identified in Ref. Q . The 
corresponding set conditions is shaded red in Fig. |3Ka). 

Particles that are much larger than the wavelength of 
light. Op ^ A, have force profiles Fz{r) that are peaked 
far enough from the optical axis that they also appear to 
be quadratic over the accessible range, but with e < 0. 
Such particles circulate in a single toroidal vortex, but in 
the sense opposite to that adopted by Rayleigh particles. 
The retrograde circulation of larger spheres was pointed 
out in a ray-optics analysis of optically trapped colloidal 
spheres [281 1 and subsequently was observed experimen- 
tally[l.‘^‘ 

Retrograde circulation in a single roll also can arise 
for particles that are intermediate in size between the 
Rayleigh range and the ray-optics regime. The example 





































force field depicted in Fig. [^[c) has this property, and is 
characterized by e < 0. Even though the quartic term is 
sizable under these conditions, the curvature of F^ir) has 
the same sign over the entire accessible range of radii, and 
only a single toroidal roll can be populated. The domain 
of such behavior is shaded blue in Fig. [3l[a). 

Topological transition and flux reversals are only pos¬ 
sible if drFz{r) changes sign in an accessible part of the 
force landscape. That occurs for conditions such as those 
in Fig.O^d), and corresponds to the green-shaded regions 
in Fig. [Sja). In such cases, the inner roll circulates in 
the same sense as the single roll in the Rayleigh regime, 
and the outer roll rotates in the opposite sense. Interest¬ 
ingly, there appears to be no set of conditions that favor 
a retrograde double-roll or the corresponding topological 
transition to a retrograde single roll. 

Different patterns of size-dependent crossovers arise for 
spheres of different materials, or in beams of different 
wavelengths or focusing properties. Figure [3] makes clear 
that small variations in properties can have a large influ¬ 
ence on Brownian vortex circulation, not simply changing 
the rate of circulation, but rather reversing the direction 
of circulation. The present work has focused on the cir¬ 
culating currents’ topology. Its results also could be used 
to address questions about drift rate and circulation fre¬ 
quency that have been raised in previous studies 0,i. 


IV. CONCLUSION 

Brownian vortexes should be generic features of 
all probability-conserving stochastic systems subject to 
time-independent driving by non-conservative forces. 
Equation (fT3]) constitutes the leading-order description 
of a weakly-driven Brownian vortex. This level of ap¬ 
proximation already captures the topological transitions 
and flux reversals that have been reported for archetypal 
Brownian vortex circulation in optically trapped colloids. 
Analytical results for topological transitions in the cur¬ 
rent density go beyond reproducing experimental results 
by clarifying their nature and providing a unified expla¬ 
nation for the system’s behavior. The idealized model of 
an optical tweezer as a harmonic well subject to a steady 
non-conservative force is likely to serve as a useful model 
other systems as well. 

It is noteworthy that the steady-state probability dis¬ 
tribution in a Brownian vortex, 

p(r)=e-'5-^W, (30) 

depends on characteristics of the force field F(r), but 
not on the particle’s mobility The particle’s transport 
properties therefore do not dictate how the probability 
redistributes as the system is driven out of equilibrium. 
In this sense, the general relationship between r(r) and 
the effective potential described in Appendix [^con¬ 
stitute an analog for the Boltzmann relation for this class 
of nonequilibrium systems. 


The generic nature of the model discussed in Sec. [TTIl 
suggests that this might be a common factor for steady- 
state circulation in a broad range of nonequilibrium sys¬ 
tems. The formalism developed here also should apply 
in other contexts such as circulatory and oscillatory flows 
in social networks, financial systems, and chemical net¬ 
works. It also would be interesting to extend this for¬ 
malism to more general systems such as systems with 
multiple fixed points, and systems of multiple interact¬ 
ing particles. 
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Appendix A: Perturbation Theory 


The effective potential, (/)(r), that determines the 
steady-state probability distribution, p(r), differs from 
the imposed potential, C/(r), by an amount 

tlj{r) = (j){r) - U{r) (Al) 

that vanishes if V x A(r) = 0. We assume therefore 
that ^/’(r) is characterized by the same small parameter, 
e, that characterizes V x A(r). Imposing conservation of 
probability through Eq. ([5|) yields a differential equation 
for '0(r), 


[VV)!-) - /3VV'(r) • VC/(r)] 

— /3 |VV'(r)|^ -f /3VV'(r) • V x A(r) 

= (A2) 

that depends on the temperature and characteristics of 

the force field, but not on the diffusing particle’s mobility. 
To first order in e, this reduces to 

V^'0(r) —• VC/(r) =/3p(r). (A3) 

The associated field 

y(r) = (A4) 


then may be obtained as an expansion in eigenfunctions 
X„(r) of the Hermitian operator 


Hir = 




(A5) 


Specifically, solutions of 

Hixn{r) = A„x„(r) (A6) 


form a complete set of basis functions labeled by index 
n with eigenvalues A„. When appropriately normalized 
they satisfy the orthogonality condition 


J Xmir)Xn{r)(fr = Smn- (A7) 
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Equation (Il3dl) follows from Eq. (IA7I) . 

Imposing conservation of probability at each order of 
the perturbation expansion ensures that solutions reflect 
steady state behavior. Higher-order corrections obtained 
by incorporating the second-order terms in Eq. (IA2I) 
will redistribute the probability distribution, but are not 
likely to eliminate qualitative features of the steady-state 
circulation that arise at first order. 


Appendix B: Generalized Lorenz-Mie theory 


has expansion coefficients, 


^mn — ^mn and 

^mn ^mnj 


(B3) 


that are related to the incident beam’s coefficients by the 
particle’s Lorenz-Mie scattering coefficients, a„ and 

For the particular case of scattering by a sphere, 
the Lorenz-Mie coefficients depend on the sphere’s radius 
and refractive index relative to the medium. 

The combined incident and scattered fields, 

E(r) =E,(r) + E,(r-rp), (B4) 


An optical tweezer can be modeled as a strongly fo¬ 
cused Gaussian beam, and its field can be expanded as a 
series in vector spherical harmonics, 

oo n 

Ei(r)=Eo^ (/sr) . 

n—1 m——n 

(Bl) 

The coefficients amn and bmn have been previously re¬ 
ported for a converging Gaussian beam [29l - IMI| . and de¬ 
pend on the numerical aperture of the lens that brings 
the beam to a focus. 

The incident beam illuminates a particle located at 
Tp, which gives rise to a scattered wave Es(r — Vp) that 
propagates to position r. The corresponding expansion 
for the scattered field, 

OO n 

E,(r)=Eo^ ^ [r™„M^^(fcr)+s™„N^^(fcr)] , 

n—1 m——n 

(B2) 


can be used to calculate the Maxwell stress tensor, T(r), 
with components 

Ttj (r) = ^mEi {r)Ej (r) - (r) Sij , (B5) 

where is the dielectric constant of the medium. The 
optically induced force is then obtained by integrating 
the Maxwell stress tensor over a surface S that encloses 
the sphere, 

F(r) h-T(r')dr', (B6) 

where n is the unit vector normal to S. Techniques to 
jj^orm this integration have been previously reported 

Figure [1] shows the computed force held acting on a 
0.5pm-radius silica sphere in water in an optical tweezer 
with vacuum wavelength A = 532 nm projected by a lens 
with numerical aperture 1.4 (convergence angle 60°). We 
dehne the particle’s point of mechanical equilibrium to be 
the origin of the coordinate system in Fig. [T] 
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